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1. Introduction 

The staggered Dirac operator, 
Jps = \Y.M^ii-K) with: (y^),^ = [/^W5,,,+^ and T7^(x) = (-l)^v<.^v^ (11) 

is the most computationally efficient way to discretize the Dirac operator, and is commonly used in 
large-scale lattice QCD simulations. Yet, it leads to 4 degenerate quark "tastes" in the continuum 
limit, and the determinant is raised to the power Nf/A in A/j-flavor simulations. The systematic 
error associated with this "rooting" is a subject of hot debate. In the literature, two different types of 
staggered-like operators which avoid rooting by representing only two tastes have been proposed: 
/) minimally doubled staggered fermions Q], which represent 2 tastes with minimal fine-tuning; 
//) staggered overlap fermions 0, where the degeneracy of the spectrum is lifted by a taste- 
dependent mass term, and the resulting operator is used as the kernel in Neuberger's overlap [||]. 

Here we study the numerical properties of (n). In spite of the additional complexity of the 
overlap, it may be simpler than the multiple fine-tuning required in (/) [^. Note also that the 
overlap kernel entering (//) can also be used without overlap, at the expense of fine-tuning. We first 
motivate the overlap kernel devised in [|T]| and study its topological properties, then consider the 
staggered overlap operator of 

2. Index of overlap kernel 

The index of a gauge field configuration is obtained from the flow of the eigenvalues of 
H{m) = 75 + rn) as a function of m. To probe the topological properties of the gluon field, it 
is essential that the varying term my^ be a taste singlet. This is not the case for staggered fermions, 
if one makes for 75 the customary choice 75 — s{x,y) = (— l)^''^''^^^, which has the decomposi- 
tion £ = 75 (g) 75 in spin (g) taste space. Instead, Adams proposes to trade e for Fj = rj^C, where 
C = ^ Li7W=penn(i234) CiCjCkCu = ^ (V^ + vj), is the symmetrized sum of 4-link parallel trans- 
porters connecting a site to its opposite in an elementary hypercube, and 'f]5{x^y) = H^^i f]^ = 
{—ly^^^^ 5x^y is the corresponding product of 4 phase factors. Fj has the spin ® taste decomposi- 
tion 75 (8) 1: it is a taste singlet, which allows mFj to probe the topology of the gauge field. 

In Figs. |l|, we compare the flow of eigenvalues with m for Adams' operator 

HA{m) = eII>s + mTi (2.1) 

and for Neuberger's operator H^/{m) = ys{Ipw +m) (where iPs and Ipw are the zero bare mass 
staggered and Wilson Dirac operators, respectively), on the same SU (3) gauge field configurations. 

The top pair of figures corresponds to the free case (the solid lines are analytic results). In the 
second set of figures, the gauge field is that of a smooth, cooled instanton. The eigenvalue flow 
shows one crossing in Neuberger's case, two nearly degenerate crossings in Adams' case, reflect- 
ing the corresponding number of flavors. The third pair of figures corresponds to a thermalized 
configuration (j(3 = 6 with Wilson gauge action). While the pattern of eigenvalue flow is less clear 
for both operators, it is remarkable that the eigenvalue gap in Adams' case seems to persist for 
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Figure 1: Flow of eigenvalues vs. m for Adams' operator (left) and Neuberger's operator (right). 



arbitrarily large values of \m\ , while in Neuberger's case additional eigenvalue crossings appear 
at negative values of m (corresponding to theories with a larger number of flavors), shifted from 
their free-field values by gauge field fluctuations. 

3. Adams' staggered overlap operator 

Having constructed a Hermitian kernel Ha (m) = elps + mT^ sensitive to the topology of the 
gauge field, one can plug this kernel into Neuberger's overlap: ^oi- = 1 + 75 sign(//(— mo)) = 
where D = y^H and = 75D75. Thus, Adams' staggered overlap operator is 

II>sov = 1 + e sign (//a (-mo)) = \ + Da/\Jd\Da, (3.1) 

with Da = eHA{—niQ) =Ips — mQeT=,. Note that the mass term —rriQeTs has the spin (g) taste decom- 
position 1<X) 75: it is "flavored". If one considers low-momentum eigenstates ^ of Ips, satisfying 

'in Adams' case, eigenvalues come in pairs X{rn) ^ — 1(— m), because eHA{m)e = —HA{—m). 
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Figure 2: Spectrum of the free Wilson operator (blue) and of the free Adams' operator (green). 
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Figure 3: Maximum magnitude maXv|Z)ot.(jc,3')| of the Dirac operator matrix elements versus Manhattan 
distance |jc — yj, for Adams' operator (left) and for Wilson's operator (right). 

Ips^ 0, they will obey (*l'^(£r5)*l') ±1 depending on their taste content. Of the initial 4 tastes, 
two combinations will give + 1 and become physical, Ught modes of the overlap operator; the other 
two combinations will give —1 and become heavy doublers. 

This is clear in the free case: Fig. § shows the spectrum of Da and of the Wilson operator Dy/, 
for mo = 1 and a free field. The splitting of the 4 tastes into 2 pairs is achieved in a symmetric way, 
more elegant than the reduction from 16 flavors to 1 in the Wilson case. Moreover, since Adams' 
kernel has a spectrum already much closer to the unit circle than Wilson's operator, one may expect 
a smaller number of operations to achieve the unitary projection D/V^D of the kernel operator 
D in Adams' case than in Neuberger's case. This is investigated in Sec. S. 



4. Locality 

First, we compare the locality of Adams' overlap operator with that of Neuberger's. In both 
cases, the matrix elements Dov{x,y) are non-zero for arbitrarily distant sites x and y, as the overlap 
operator is not ultra-local. What matters, however, is the decrease in magnitude of |Dov(.'c,y)| with 
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Figure 4: Computer cost of one overlap propagator, measured in outer CG iterations (left), matrix-vector 
multiplications (middle) and CPU time (right). The values for Adams' operator are shown in red, those for 
Neuberger's operator in green. The gauge field is the free field (top), or a j3 = 6.0 configuration (bottom). 



the distance — 3^1, which should be bounded by exp(— |x — j|/(ca)), where {ca) is a localization 
length proportional to the lattice spacing a, and thus shrinking to zero in the continuum limit. 

Fig.^shows the maximum magnitude max^, |Dov(jco,3')| versus the Manhattan distance \xq —y\, 
chosen to follow the conventions of Ref . for the Neuberger operator. The left figure corresponds 
to Adams' overlap operator, the right one to Neuberger's, on the same gauge configurations at 3 
values of jS. While Adams' operator behaves differently at short distance because of the 4-Iink 
transporters, at large distances the decay of the matrix elements is exponential as in Neuberger's 
case, with a similar localization length. ^ 

5. Robustness to gauge fluctuations and efficiency 

We have studied the computational cost of a quark propagator calculation with Adams' op- 
erator and compared it to Neuberger's propagator (for one component) on the same gauge field 
background, and with the same numerical approach. As the matrix to invert is 4 times smaller, and 
its spectrum is closer to the unit circle, at least in the free case (Fig. Adams' operator may be 
computationally cheaper. 

The propagator is obtained as the solution of {Dm, + my{Dov + m)x = {Dgy + m)^b, using a 
conjugate gradient iterative solver, using the following simple and robust method [^: at each 
iteration of this outer CG, sign(//) is applied to a vector v through a Lanczos process, building a 

^This happens even though the kernel of Adams' operator is much less local than that of Neuberger's: a less 
ultralocal kernel may lead to a more local overlap operator [0]. 
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Figure 5: Spectra of Adams' operator. Left to right: free field, j3 = 6.0, j3 = 6.0 with larger p, j3 = 5.8. 



tridiagonal matrix. Its eigenvalues are representative of those of H, and we replace them by their 
sign. The results are presented in Fig. for a free field (top) and a j8 = 6 configuration (bottom). 
The 3 figures in each row show the relative norm of the residual, |r|/|ro|, v^'. outer CG iterations 
(left), number of matrix-vector multiplications (middle) and CPU time (right). 

In the free field case, the CPU time to find the solution for Adams' operator is almost an order 
of magnitude smaller than for Neuberger's, thanks to a construction of the sign function requiring 
fewer matrix-vector multiplications, each with a smaller CPU cost. The inversion of the unitary 
operator converges at the same rate, reflecting the similar spectral properties in the infrared. 

The situation changes on a j8 = 6 configuration. The outer CG now converges noticeably faster 
in Neuberger's case (left). This advantage is offset by the cost of the sign function, which still is 
cheaper in Adams' case (middle). Finally, the CPU time per matrix-vector multiplication is a factor 
ff{2) smaller in Adams' case (Adam's matrix is one quarter the size of Neuberger's, with each site 
connected to (8-1-16) neighbours, instead of 8 with 2 Dirac components). In total, the CPU time to 
find the solution is only &{2) times smaller in Adams' case. 

This loss of efficiency can be traced to changes in the spectrum of Adams' kernel Da in the 
presence of gauge field fluctuations, as illustrated in Fig. ||. While the free spectrum (left) is 
remarkably close to the unit circle, it is quite different at j3 = 6 (2nd panel). The splitting of the 4 
tastes into 2 pairs is markedly reduced. The reason is that the taste-dependence of the mass operator 
ePs is achieved via 4-link transporters: fluctuations in the gauge links are raised to the 4th power. 
It is the same reason for which the chirality (*F^r5*F) of near-zero modes of the ordinary staggered 



operator is so small []10|]. Here, one may attempt to restore the mass splitting of the pairs of tastes, 
by increasing the mass parameter mo, called p in Adams' Ref. [Q]. The effect of such increase is 
shown Fig. |5| (3rd panel). The gap in the spectrum, which was the complete unit disk in the free 
case, and which shrank to a small but disk-like shape at jS = 6 for wiq = 1, now becomes a very 
narrow band. Unitary projection of the operator becomes more difficult, and after unitarization 
many modes are present near the origin, which makes inversion more difficult as well. 
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This figure also shows that the spectrum remains centered about the origin: changing mo is 
not the analogue of changing the mass in the Wilson operator, which shifts the whole spectrum. 
Rather, mo is the analogue of Wilson's hopping parameter r 

This is why the eigenvalue gap in the Hermitian operator HA{m) (Fig. |l]) persists for large val- 
ues of |m| . Shifting the whole spectrum of Da by a taste-independent mass term is also possible, but 
will destroy the symmetry of the spectrum about the origin without any computational advantage. 
Finally, Fig. || (right) shows how the gap in the spectrum of Da closes at j8 = 5.8. 

6. Conclusion 

Our study shows that Adams' staggered overlap operator works as advertised: the taste- 
dependent mass operator in its kernel yields 2 massless tastes without fine-tuning, and the topology 
and locality properties are similar to Neuberger's operator. 

On very smooth gauge configurations, the computer cost of a quark propagator is nearly an 
order of magnitude less than in Neuberger's case, but the 4-link transporters in the flavored mass 
term reduce this advantage to a factor 0'{2) on p = 6 configurations. Another drawback of Adams' 
construction is that the continuous symmetry of the massless staggered overlap operator is U{1), 
not SU (2) as one would wish for a 2-flavor chiral symmetry. 

The lack of robustness and of full chiral symmetry can both be addressed by modifying the 



mass operator, for example with 2-link transporters [11 1 that reduce the number of light tastes to 



1. Preliminary results [|12|| confirm our expectations, but do not bring the cost of staggered overlap 



fermions near that of ordinary staggered fermions: avoiding rooting still has its price. 
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